Abstract 



We study the properties of the probability density function (PDF) 
of a bistable system driven by heavy tailed white symmetric Levy 
noise. The shape of the stationary PDF is found analytically for the 
particular case of the Levy index a = 1 (Cauchy noise). For an arbi- 
trary Levy index we employ numerical methods based on the solution 
of the stochastic Langevin equation and space fractional kinetic equa- 
tion. In contrast with the bistable system driven by Gaussian noise, 
in the Levy case the positions of maxima of the stationary PDF do 
not coincide with the positions of minima of the bistable potential. 
We provide a detailed study of the distance between the maxima and 
the minima as a function of the potential's depth and Levy noise pa- 
rameters. 
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1 Introduction 

The term "Levy noise" stands for a class of non-Gaussian random 
stationary processes possessing alpha-stable Levy probability distri- 
butions, the class of probability laws originally studied by French 
mathematician Paul Pierre Levy [I]. These laws are ubiquitous in 
nature due to generalized central limit theorem [2]. The alpha-stable 
probability density functions (PDFs) exhibit slowly decaying power- 
law asymptotic behavior of the form | a? | 1 ^ , where a is called the 
Levy index, < a < 2; therefore, the variance is infinite. For that 
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reason Levy noises naturally arise in the description of random pro- 
cesses with large outliers, far from equilibrium, see, e.g., [21 IU G3 E] 
and references therein. 

The study of relaxation processes in dynamical systems sub- 
jected to the Levy noise is of interest for several reasons. From the 
nonequilibrium statistical physics point of view such systems represent 
paradigmatic examples of non-Brownian processes in the systems far 
from equilibrium. Also, it has been demonstrated that such systems 
can be used as "minimal" models for the description of such diverse 
phenomena as anomalous transport in turbulent plasmas [7, 8, 91 HOj 
and abrupt climatic changes \12\ [T3]. Another important field of 
research is related to the reverse engineering problem |14| . 

The properties of linear dynamical systems subjected to the Levy 
noise have been studied in the two seminal papers |15[ [TC] and later 
in |17l IT8] . The relaxation in nonlinear systems has been extensively 
studied within the last decade, see the reviews [5j[6l[T9] and references 
therein, and more recently, e.g., in [201 l2T| l22j l23l 124] [25] . 

In the present short communication we demonstrate the first re- 
sults on stationary states in nonlinear dynamical system with a sym- 
metric double-well potential. We find analytical expression for the 
stationary PDF for the Cauchy case a = 1 and employ numerical 
methods to investigate properties of the PDFs for arbitrary a. 



Our starting point is the stochastic Langevin equation which can be 
written in dimensionless variables as 



where U(x) = x 4 /4 - ax 2 /2, £ a (t) is an a-stable noise source with 
the unit intensity, a is the noise amplitude, a > is the parameter 
that allows to govern the potential's well depth (note, that ±\fa are 
the positions of its minima). With such a choice of dimensionless 
variables we can compare our results with those obtained for a quartic 
potential |26j, setting a = 0. This kind of a potential, along with the 
Langevin approach framework is used to describe a large variety of 
systems where a switching between two metastable states occurs, see, 
e.g. Refs. [HUES 03] mentioned above. The corresponding to Eq. ([!]) 
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Fokker-Planck equation reads as [T71 [TH] 



df _ d fdU \ ff»f 

dt dx \ dx J ) d\x\ a ' y ' 

where d a f /d\x\ a is the Riesz fractional derivative |27j . 

We consider the stationary case. For the characteristic function 
defined as f(k,t) = f^° 00 dxf(x,t)e lkx we get: 

g + a f-„W-/ = 0. (3, 

If we restrict ourselves to the case a = 1, the latter equation simplifies 
to the following: 

Looking for the solution in the form f(k) = Ce zk , we ar- 
rive at the characteristic equation z 3 + az — a = 0, for which 
the three roots are determined by the Cardano formulas: Z\ = 
-(p + q)/2 + i\/3(p - q)/2, z 2 = z*, Z3 = p + q, where p = 

a/2 + ^cj 2 /4 + a 3 /27) 1/3 , q = (a/2 - ^/a 2 /4 + a 3 /27^ ^ , and 

the general solution is written as f(k) = C\e zk + C^e 2 k + C^e Z3k , z = 

The coefficient C3 is zero due to f(k) — > 0, k — > 00. The coef- 
ficients C\ and C2 are determined from the normalization condition 
f(k = 0) = 1 and the boundary condition df(0)/dk = arising due to 
symmetry of the problem. Thus, for the stationary state we have 

f(k) = j^(ze**\ k \-z*e*W). (5) 

Making an inverse Fourier transform we get the stationary PDF in 
case a = 1 : 

/(X) = ~~Vx 4 + {z 2 + z* 2 )x 2 + |z| 4 ' (6) 

The minima of the U(x) are located at ixjj, where xjj = 
yfa. The maxima of the stationary PDF are located at ±Xf, 
where Xf = ^—(z 2 + z* 2 )/2 (these maxima always exist due to 

z 2 + z * 2 = - \(^4 a y 2 7 + a 2 -a) 2/3 + (^4^/27 + a 2 + a) ^ + 
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2 8 / 3 a/3 /2 2 / 3 < 0). Clearly, x f and xjj are not the same. This fact 

was expected from the previous studies [26], where the bimodality of 
the stationary PDF for quartic potential was revealed. 

Let us introduce a parameter A = (xf — xjj)jxu. In the weak 
noise limit, a 3 /a 2 3> 1, we get A oc a 2 /a 3 . In the strong noise limit, 
a 3 /a 2 <C 1, one may become convinced that Xf ~ er 1 / 3 /y/2 and A ~ 

1 /6 

(cr 2 /a 3 ) /v2 — 1. The last expression reveals that the distance 
between the corresponding minima and maxima becomes unlimited 
with the increasing noise intensity. 



3 Numerical results 

In this section we employ two different techniques to simulate the 
stationary PDFs. The first method is based on the numerical solution 
of space- fractional Fokker-Planck equation by means of the Griinwald- 
Letnikov approximation of the Riesz fractional derivative [28]. The 
second method uses direct Langevin dynamics simulations. 

Since both these approaches give us stationary PDF it is straight- 
forward to compare their results to each other, as well as to the an- 
alytical results obtained in Section [2] In the kinetic approach we 
investigate the system by varying potential's parameter a, while fo- 
cusing on the noise intensity properties in the Langevin simulations 
instead. 

As the first test of the approaches, let us compare the analytical 
results obtained in Section [2] with the corresponding numerical data 
of the Langevin simulations as it is shown in Fig. [T^,. Evidently, the 
simulations (dots) do match well with the analytics (solid lines) at 
all noise intensities. The Figs. [TJd and [Tj: show the results of both 
simulation methods, the Langevin (points) and the kinetic (solid or 
dashed lines) for different values of the Levy index. 

The double logarithmic scale in Fig. [TJd makes it obvious that the 
slopes of the asymptotes have a disrepancy which is probably due to 
a growing inaccuracy in extracting PDF from of the Langevin simula- 
tions at long distances. Fig. [Tp gives a quantitative comparison of the 
methods' results plotting the stationary PDF peaks f{xA and min- 
ima /(0), and the asymptotes slopes slope = (—din //din x), x — > oo 
showing a good agreement between the results obtained with the use 
of both simulation techniques. In Fig. [TJ; we also plot two analyti- 
cal results. The values of f(xj) and /(0) computed from Eq. (roj) are 
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Figure 1: (colour online) (a): Comparison of analytical (solid line, Eq. [6]) 
and numerical (dots) results for stationary PDF in double-well potential, 
a = 1, a = 1. (b): Comparison of Langevin and kinetic simulations results. 
Shapes of PDFs in linear and log- log scales, a = a = 1. (c): Quantitative 
comparison of the shape parameters, a — a — 1. 

shown by two black crosses. The blue solid line shows the asymptotic 
slope of the stationary PDF f(x) oc l/|x| a+3 (see Eq. (71) in [2H]) ob- 
tained for a quartic potential well. The discrepancy between this slope 
and those obtained with numerical simulations (dashed line and trian- 
gles) indicates the necessity of taking into account the quadratic term 
while estimating the asymptotic properties of the stationary PDF in 
the double-well potential. 

In Fig.[2]the shapes of stationary PDFs are shown for different val- 
ues of the parameters a, a and a. The main feature of the stationary 
PDFs is demonstrated here clearly: the PDF's peaks do not coincide 
with potential's minima. The results of quantitative analysis of this 
phenomenon is presented in Fig. [3] in terms of the coefficient A. As 
it should be, A — » at a — ^ 2, however, A ^ at other Levy indices 
of the noise. The deviation becomes more visible at higher noise in- 
tensities and when the potential's minimum becomes more shallow. 
Interestingly, our numerical simulations demonstrate the existence of 
maximum at intermediate a's in Fig. [3^,. Both figures show a good 
agreement of the results obtained numerically and analytically. 

4 Conclusions 

Our studies of the Levy flights in a double-well potential demonstrate 
that the positions of the maxima of the stationary probability density 
functions do not coincide with the positions of the minima of the 
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Figure 2: (colour online) (a,b): Langevin simulations of stationary PDFs 
in linear and log-log scales, respectively. Case a = 0.5, a = 1. (c,d): Kinetic 
simulations of stationary PDFs. (c): a = 1.05, a = 1.0, a = 1; (d): a = 1.50, 
cr = l. Arrows show the positions of the potential minimum's location yfa 
(1) and the PDF's peak coordinate Xf (f). 
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Figure 3: (colour online) Quantitative analysis of the distance between the 
PDF peaks and potential minima, (a): The results of Langevin simulations 
for various noise intensities a and a = 1 are shown with filled squares and 
triangles; with empty □ we depict the data obtained via the numerical kinetic 
approach for a = 1. (b): The results of the kinetic approach for various 
coefficients a alongside with analytical results marked with x . 
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potential well. In the present short communication we provide a first 
study of the properties of the stationary PDFs and reveal that the 
difference between the mentioned maxima and minima grows with 
the increasing noise intensity and decreases with the potential's depth. 
Very recently, the discrepancy issue between the shape of the potential 
and the shape of stationary PDF was addressed from the point of 
view of "mean reversion" and "mode reversion" [29] . Interestingly, the 
same discrepancy was found there under certain conditions even in the 
Gaussian case. We plan to elaborate on the "reversion" phenomena 
for Levy flights in the forthcoming long paper. 
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